Quasiparticle thermal conductivity in the vortex state of high-T c cuprates. 
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We present the results of a microscopic calculation of the longitudinal thermal conductivity, k xx , of a 
d-wave superconductor in the mixed state. Our results show an increase in the thermal conductivity 
with the applied field at low temperatures, and a decrease followed by a nearly field independent 
k xx (H) at higher temperatures, in qualitative agreement with the experimental results. We discuss 
the relationship between the slope of the superconducting gap and the plateau in k xx (H). 



o 
o _ 

u • 
Of 

s : 



i 

-a 
c 

o 
o 



(N 
> 

O 

o 

o\ 

-i— > 



C 

o 
o 



X 



While the origin of the unusual normal state properties 
of the high-T c cuprates remains controversial, many ex- 
perimental results lend support to the belief that the su- 
perconducting state of these materials is relatively mun- 
dane, and consists of a pair condensate with an "uncon- 
ventional" order parameter, believed to be of d-wave sym- 
metry, and well-defined quasiparticle excitations above 
it. However, transport measurements in the vortex state 
yield some surprising results which question whether the 
low-tempcrature physics of the cuprates is well under- 
stood. 

Perhaps the most notable such example is the observa- 
tion of a plateau-like feature in the magnetic field depen- 
dence of the longitudinal thermal conductivity, k xx (H), 
in superconducting Bi 2 Sr 2 CaCu208+5 (BSCCO) Q and 
YBa2Cu307_,5 (YBCO) Q. At temperatures above a few 
Kelvin k xx (H) initially decreases with applied magnetic 
field, H, before becoming nearly field independent. Q-^] 
The kink and hysteresis in n xx (H) observed in zero-field 
cooled BSCCO @|§ have not been found in YBCO |] 
or in the field-cooled BSCCO samples Q], and we take 
the point of view that they are not a generic feature of 
the k xx (H) dependence in cuprates. In some samples 
the "plateau" is not reached, and k xx (H) decreases con- 
tinuously for H < 14 Tcsla, ||,|| as it does for conven- 
tional superconductors for H <C H C 2- || This is in sharp 
contrast to the behavior of the thermal conductivity of 
the cuprates at ultra-low temperatures, which increases 
with H. Moreover, both features have been observed 
in the same sample ||, suggesting that the two regimes 
should be understood within the same approach, which 
presents a challenge to theories of quasiparticle transport 
in the vortex state of the cuprates. 

Existing theories are based on the observation that, for 
a single vortex, the properties of superconductors with 
lines of nodes in the energy gap are dominated by ex- 
tended quasiparticle states, rather than, as in the s-wave 
case, by localized states in the vortex cores. The mag- 
netic field is taken into account semiclassically, by in- 
troducing a Doppler energy shift due to the circulating 
supercurrents. In one approach physical quantities are 
computed locally and are averaged over a unit cell of the 



vortex lattice [|3| . The results obtained within this frame- 
work describe experimental results quite well |q-^0[; and 
the predictions for k xx (H) by Kiibert and Hirschfeld || 
fit the ultralow-T measurements. |(| The results of Ref. 
H show a decrease (H) at higher T <C T c . but 

not a plateau. A similar approach was used for heavy 
fermion superconductors. |]llf 

Franz [[l2| has pointed out that the averaging pro- 
cedure used in Refs. |§[-|lC|] neglected scattering from 
the vortex lattice. He included the scattering semi- 
phenomenologically, assuming a completely disordered 
vortex lattice and a Matthiesen type rule for the impu- 
rity and vortex contributions to the lifetime |13|. He pre- 
dicted that at any T <T C k xx (H) will approach asymp- 
totically the universal value 0Jl^l koo = vrTA u 2 /6A , 
where Nq is the normal state density of states at the 
Fermi surface, v is the Fermi velocity, and Ao is the gap 
amplitude. His results are in qualitative agreement with 
the results of Ref. || , but the theory requires additional 
assumptions to explain the low-T limit. Consequently, 
no existing theory explains both the low-T results and 
the plateau. 

In this work we present a microscopic theory of the 
thermal conductivity of a two-dimensional clean (I ^S> £o, 
where I is the mean free path, and £o is the coherence 
length) d-wsive superconductor in the vortex state, with 
the field perpendicular to the basal plane, and thermal 
current in the plane. Our results fit the low-T data of 
Ref. H and show a plateau-like feature at higher T. We 
employ the quasiclassical method jl6|,[l7|], in which the 
basic quantity is the matrix Green's function integrated 
over the quasiparticle energy, 



-^G(p,R;w„,uv) = ifi 
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Here the diagonal elements, g and g, are the particle 
and hole propagators, respectively, the off-diagonal ele- 
ments, Gorkov's anomalous functions / and f>, are prob- 
ability amplitudes for the creation or destruction of a 
Cooper pair, R is the center of mass coordinate, £ p is 
the energy of a quasiparticle with momentum p mea- 
sured relative to the Fermi surface, and 4> is the angle 
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parameterizing the cylindrical Fermi surface. The spa- 
tial dependence of the vortex lattice (VL) is modeled by 
an Abrikosov-like solution A(R, 0) = A(R)ip(4>), where 
A (R) = Efe H C ky e lk yyeM-(x- A 2 k y ) 2 /2A 2 ), the mag- 
netic length A = (2eH)^ 1 ^ 2 ( h = c = 1 throughout the 
paper), and <p(<ft) — cos 20 is the angular dependence of 
the gap. g§ 

We use the quasiclassical method in the framework 
proposed in Ref. ]2C| ] for s-wave superconductors near 
H C 2- In that approach the diagonal elements of g are 
approximated by their spatial averages, while the exact 
spatial dependence of the off-diagonal elements is kept. 
The authors of Ref. ^(J noticed that g and g are periodic 
in R with the periodicity of the vortex lattice, and that 
Fourier components (?k with reciprocal vortex lattice vec- 
tors K 7^ are smaller than 3k=o (which is the spatial 
average of g) by a factor exp(— A 2 i^ 2 ). This approach 
provides a very accurate description of the vortex state 
in conventional superconductors. 20]-24 The analysis of 
Ref. [^O) relies on the fact that the spatial dependence 
of the order parameter is given by an Abrikosov-like so- 
lution of the linearized Ginzburg-Landau equations, and 
therefore works especially well in extreme type-II super- 
conductors where it can be extended over almost the en- 
tire region of linear magnetization. |25j ] For an s-wave 
superconductor it breaks down at fields where proper- 
ties are determined by the states bound to vortex cores. 
In unconventional superconductors where properties are 
dominated by extended states the approach works 
down to fields close to H c i> i-e. for all experimentally 
relevant fields. jl9| 

The leading order quasiclassical propagator is deter- 
mined from equations [^3| , p6| 

(2w„ + v(V - 2ieA))f = 2iA(K)g cos 2<f> (1) 
(2S„ - v(V + 2ieA))P = 2iA*(R)g cos 20, (2) 

complemented by the normalization condition g 2 — // T = 
— 1. Here A = (0,Hx,0), and the renormalized fre- 
quency, iui n — iu> n + ia(iuj n ), depends on the self-energy 
due to impurity scattering, er, which is determined within 
the t-matrix approximation p7[ . We ignore corrections 
to the impurity vertex due to the magnetic field, which 
is justified when A <C I. Since A ~ 183A/VA, Tesla, for 
clean samples this condition is satisfied at least above 
0.1-0.5 T. 

The solution of Eqs.(Q),(|) is known to be [||j2C| 

g = -isgn(w„)P(0,w„), where P(0,w„) = [1 - 
i^Tr(2AA/v) 2 W'(u n )cos 2 20]" 1 / 2 , the function W{u) = 

e~ u erfc(— iu), u n = 2za3 ra Asgn(cj„)/w, and A is the spa- 
tial average of the amplitude of the order parameter. In 
the limit H = setting A equal to the gap amplitude, Ao, 
we recover the BCS Green's function and consequently 
all the results of the standard "dirty d-wave" theory. p§| ] 
For a pure sample at low fields the residual density of 
states (DOS) N(0)/N o = (u/AA 7r\/2) ln(8v/2AA /u) oc 



V H \i\{H C 2 1 H); the field dependence differs by a loga- 
rithm from the single vortex result of Volovik 0. For 
the range of fields where comparison with experiment 
has been made the logarithm is nearly constant, and the 
result for the DOS agrees with that of of Ref. ^ up to a 
numerical factor. In the following we consider scattering 
in the unitarity and Born limits, so that the parame- 
ters of the theory are the normal state scattering rate, 
F, and the dimensionless quantity AA/v. The latter can 
be estimated by setting A ~ A (l — H/H^) 1 ^ 2 , and 
noticing that AAq/v = (Akf /2)(v2/v), where V2 is the 
slope of the gap at the node. In BSCCO kf ~ 0.737 
A -1 , and is relatively doping independent, |29[] then 
AA/v ~ cty/(l - H/H c2 )/H; where, for field in Tesla, 
a ~ 3.5 for the overdoped case, where vlvj. ~ 20, but 
becomes smaller on the underdoped side. (30j In YBCO, 
taking k f ~ 0.8A" 1 |n|) we obtain a ~ 5 if v/v 2 ~ 14 
||, or twice that value if v/v2 ~ 7. [{32[ 

The heat current can now be determined from the lin- 
ear response equations (24|]2(| 

j h =«N lim V /"#v(0)(v(0) ■ VT) (3) 

„ i ^ g(u n ) - g{u n - t^o) 

xw„ w„ - uj ) — — — — , 

lUJo + 26 

where b = (Ap(u> n ) + A*f(oj n - lo )) / (g(uj n ) - g(uo n - 
u>o))- This leads to the thermal conductivity 
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Since Re[P(u, 4>)] is the angular dependent DOS, and 
t(2, tp) is a transport lifetime, Eqs.([|) and ^ generalize 
the normal state expression k xx = K n = tt 2 Nqv 2 Tt/6. 
There is a contribution to the quasiparticle lifetime from 
Andreev scattering off the vortex lattice, which has the 
same symmetry on the Fermi surface as the gap function, 
and is therefore the more important the higher the quasi- 
particle energy. This contribution appears even for a pe- 
riodic vortex lattice, since Bloch's theorem only forbids 
scattering from a periodic potential in the absence of im- 
purity broadening of the quasiparticle wave-packet, not 
in the presence of impurities. Note that the dependence 
of the DOS and the scattering rate on angle cannot be 
disentangled, in contrast to the assumption of Ref. [12| . 
Also, it is clear from the form of Eq.(^) that a simple 
scaling form as a function of T / \fll cannot be expected. 

First, we consider the limit T = 0, where analytic re- 
sults can be obtained. It is convenient to analyze the 
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thermal conductivity in terms of a transport lifetime th, 
defining k xx {H)/T = ttNqv tjj /6, and expressing tjj in 
terms of the lifetime r e // = 1/(27), where 7 = cr{u> = 0) 
is the zero-frequency scattering rate. For H = we ob- 
tain tjj = kE(k)/TrAo, where E is the complete elliptic 
integral, and k = (1 + (2Aor e //)~ 2 ) -1 / 2 in agreement 
with Refs. ||^,[l5|]. For Aor e // > 1 we recover the uni- 
versal limit k xx = Koo 

= t:TN v 2 /6A . At finite 

fields the integral in Eq.(Q) is evaluated approximately 
to give tjj ~ (T e //i;/4 % /7rAA§) 1 / 2 . In the "clean" limit 
v/A ^> 7 for the important cases of Born (r^) and uni- 
tarity (IV) scattering, we find 

_ f (£)V4 (rsAo) -l/2 ln -l/2 ( 8^ ); 

TH \ (^-^(IVAo)- 1 / 2 ^) la 1 ^!^). 

(6) 

This expression is in a good agreement with the numeri- 
cal results at intermediate fields. In Fig. 1 we use the nu- 
merical evaluation of Eq. (||) to compare k xx (H) for Born 
and unitarity scattering. In the case of Born scattering 
k xx {H) only reaches the universal limit at exponentially 
small fields, and is approximately field-independent at 
higher H, Eq.(^). In the unitarity limit, on the other 
hand, k xx (H) increases approximately as \/~H at low 
fields and the approach to k 00 is clearly seen. At low 
T it is seen that k xx (H) tends to koo over a field range 
of order a few Tesla ||, and that the field dependence 
is close to \f~H at low T. Q This suggests that impu- 
rity scattering is close to the unitarity limit at low tem- 
peratures, and we consider this limit hereafter. The in- 
dependence of tjj differs from that found in the Doppler 
shift theory Jj|, however, since varies slowly, for 

experimentally accessible fields the behavior of k xx (H) 
is close in to that obtained by Kiibert and Hirschfeld. ||] 
In Fig. 2 we fit the data of Ref . |(| . Using the the value 
of v/v 2 = 14 from Refs. we find V to be about 4 

times greater than the estimate of Ref. || . In the present 
theory k xx {H) increases faster with field the smaller the 
slope of the gap, approximately as yjH/a 2 Y, see Eq. (||), 
so that the same data can be fit on assuming a smaller 
impurity concentration if we allow for a larger ratio V2/V. 
|p2| This also implies that a smaller slope of the gap can 
result in a larger increase in k xx (H) for a dirtier sample, 
see Fig. 2, which may be a part of the explanation of 
a faster increase of k xx {H) in nominally dirtier BSCCO 
compared to YBCO. § 

At finite temperatures the integration in Eq.([I|) has to 
be carried out numerically; we confine ourselves to the 
regime where the gap is temperature independent. Wc 
have checked and find agreement with previous results 
at H — 0. |lj] k(0)/T increases rapidly for T ^ 7, and 
for clean samples at T ~ 0.1 Ao it exceeds koo by two 
orders of magnitude. At the same time at strong fields, 
»/A » T, k xx (H)/T is controlled by the field strength 



and is only weakly temperature dependent. The field de- 
pendence at T — is weak, close to VH, so that for 
H ~ 10 T k xx (H)/T < IOkoo/T. An immediate conclu- 
sion is that for T > 7 there is a steep drop in k xx (H) 
at low fields, followed by a crossover to the \/~H like be- 
havior. This crossover occurs when v/A ~ T, at a field 
H* oc T 2 , as indicated experimentally. p]|| 

In Fig.3 we show the results: at high T the ratio 
Kxx(H) / k xx (0) decreases continuously as a function of 
field, however as the temperature is lowered a plateau 
develops, and at even lower T a pronounced minimum ap- 
pears in the field dependence. As the lower panel demon- 
strates there is a transition to an increasing k xx (H) at 
the lowest temperatures. Numerically, we find that at 
fixed T a smaller value of a corresponds to a more pro- 
nounced plateau (H). This implies that the nearly 
field independent thermal conductivity should be more 
evident in the underdoped cuprates. []30"f 

To summarize, we give a theory of the thermal con- 
ductivity of a (i-wave superconductor in the vortex state 
as a function of field at T -c T c . In zero field and at low 
T we recover the universal results |l4|,[l5|]. We find that 
there exists a contribution to the transport lifetime due to 
scattering off the vortices, which has the same symmetry 
as the gap function. Consequently, at low temperatures 
when most of the quasiparticles are near the gap nodes it 
plays a minor role and the effect of the field is to increase 
the DOS and the thermal conductivity. In contrast, at 
high T the main effect of the vortices is to introduce a 
new scattering mechanism, and n xx {H) decreases sharply 
to a plateau-like feature at higher fields (a broad mini- 
mum); the crossover field scales with T 2 . This behavior 
is in agreement with experimental results. Fi- 
nally, it is important to compare the predictions of our 
theory with the measurements of the transverse (Hall) 
thermal conductivity, n xy , 0] to which phonons do not 
contribute. At low fields n xy is given by Eq.(Q) with the 
kernel T xy (oj, 4>) = uj c t(uj, cj))/2Re[a(uj)]. However, as the 
measurements are taken at T > 20K, in a regime where it 
is necessary to include inelastic scattering in the analysis 
of both longitudinal and transverse thermal conductivity, 
we defer a discussion of this issue until later. 
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FIG. 1. Field dependence of the thermal conductivity at 
T = for scattering in the unitarity (solid lines) and Born 
(dashed lines) limits. Inset: low field behavior, normalized by 
the universal value kqo- H C 2 = 150 T. 
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FIG. 2. Field dependence of k„ (H ) at T = compared to 
data from Ref. [6]. 
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FIG. 3. Top: field dependence of k X x(H)/k(0) for 
r = 0.003Ao and a = 7 at temperatures where n xx (H) is 
decreasing; bottom: the change in the thermal conductivity 
as a function of field for different temperatures. 
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